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We report on the first direct search for charged Higgs bosons decaying into cs in tt events produced 
by pp collisions at y/s — 1.96 TeV. The search uses a data sample corresponding to an integrated 
luminosity of 2.2 fb _1 collected by the CDF II detector at Fermilab and looks for a resonance in the 
invariant mass distribution of two jets in the lepton+jets sample of tt candidates. We observe no 
evidence of charged Higgs bosons in top quark decays. Hence, 95% upper limits on the top quark 
decay branching ratio are placed at B(t — > H + b) < 0.1 to 0.3 for charged Higgs boson masses of 
60 to 150 GeV/c 2 assuming B(H + — > cs) = 1.0. The upper limits on B{t — > H + b) are also used 
as model-independent limits on the decay branching ratio of top quarks to generic scalar charged 
bosons beyond the standard model. 
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The standard model (SM) is remarkably successful in 
describing the fundamental particles and their interac- 
tions. Nevertheless, it is an incomplete theory. An im- 
portant unresolved question is the mechanism of elec- 
troweak symmetry breaking (EWSB). In the SM, a sin- 
gle complex scalar doublet field breaks the symmetry, re- 
sulting in massive electroweak gauge bosons and a single 
observable Higgs boson To date, the Higgs boson has 
not been discovered, and consequently the mechanism of 
EWSB remains in question. 

Beyond the SM, many diverse hypotheses with ex- 
tended Higgs sectors have been proposed to explain 
EWSB. The simplest extension is a two Higgs-doublet 
model (2HDM). The minimal supersymmetric standard 
model (MSSM) employs the type-II 2HDM, where at 
leading order one doublet couples to the up-type fermions 
and the other couples to the down- type fermions 0. 
The two Higgs doublet fields manifest themselves as 
two charged Higgs bosons (H^) and three neutral Higgs 
bosons (h, H, A). 

In 2HDM and MSSM, the top quark is allowed to decay 
into a charged Higgs boson (H + ) [3J] and a bottom quark. 
The tree level branching ratio of top quarks to H + , B(t — > 
H + b), is a function of the H + mass {m H +) and tan/3. 
The parameter tan (3 is the ratio of vacuum expectation 
values of the two Higgs doublets. In MSSM, B(t ->■ H+b) 
also depends on extra parameters related to the masses 
and couplings of the other supersymmetric particles. The 
B(t — » H + b) is relatively large if tan/3 is small (< 1) or 
large (> 15) At low tan/3, H+ predominantly decays 
into cs for low m H + (< 130 GeV/c 2 ) and t*b (-*• Wbb) 
Q for higher m H + . In the high tan (3 region, the H+ 
decays into t+v almost 100% of the time. 

At Tevatron collider experiments, H+ searches have 
been performed for the H+ — > tv in ti decays. Some 
searches placed direct upper limits on B{t — ► H + b) by 
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taking advantage of the expectation that B(H + — > tv) 
— 1.0 at high tan/3 [(|. Other searches set limits on the 
MSSM parameter plane (m H + , tan fi) using inclusive H + 
decay branching ratios in the MSSM Q. The various H+ 
final states supplement the SM ti decay channels. The 
previous searches focused on measuring deviations from 
the SM prediction for the ti production and decay, rather 
than reconstructing H+ bosons. 

In this Letter, we report on the first direct search for 
H+ — > cs produced in top quark decays by fully recon- 
structing the cs mass. The final state of H+ — > cs is 
mostly two jets, as is the hadronic decay of the W boson 
Q in SM top quark decays. The search is performed by 
looking for a second peak in the dijet mass spectrum (in 
addition to that from the W boson) in top quark decays. 

In the SM, each top quark decays into a W boson and 
a 6-quark exclusively. In this analysis we use the lep- 
ton+jets ti sample [9(, where in the SM one W decays 
to quarks (qq 1 ) and the other W decays to ev or [iv. 
Each final-state quark is assumed to form a hadronic jet; 
the jets are clustered using a cone algorithm with a cone 
radius AR (= ^(Arj) 2 + (A0) 2 ) of 0.4 [10]. This lep- 
ton+jets sample has a good signal-to-background ratio 
for tt and is ideal for dijet mass analysis. 

The CDF II experiment at the Fermilab Tevatron mea- 
sures the products of proton-antiproton collisions at ^/s 
= 1.96 TeV. The lepton momentum is measured using 
an eight-layer silicon microstrip detector and a cylindri- 
cal drift chamber immersed in a 1.4 T magnetic field. 
The energies of electrons and jets are measured using 
calorimeters with acceptance up to pseudorapidity as of 
\r/\ = 3.6 11]. Charged particle detectors outside the 
calorimeter identify muon candidates up to \rj\ = 1.0. 
Details of CDF II can be found elsewhere [12]. 

Lepton+jets ti events are selected by requiring an elec- 
tron or a muon with px > 20 GeV within |?7|=1 and by 
requiring missing transverse energy larger than 20 GeV 
to account for the neutrino [HI]. Then, the four most 
energetic jets (called leading jets) within \r/\ < 2.0 are 
required to have Et > 20 GeV after jet energy correc- 
tions [10]. In addition, at least two of the leading jets 
are required to contain a long-lived hadron containing a 
6-quark [3] by demanding that these jets contain tracks 
forming a displaced secondary vertex (called a 6-tag). 

The SM processes are regarded as backgrounds for the 
H + search. The largest background is W bosons in SM ti 
events (92% of the total background) . The rest of the SM 
processes are referred to as non-tf backgrounds. These 
include W+jets, multijets, Z+jets, diboson (WW, WZ, 
ZZ). and single top events. Details of the non-tt back- 
ground estimation method are given in [3] ■ Assuming a 
ti cross section of 6.7 pb [IH and a top quark mass of 175 
GeV/c 2 , we expect 152.6 ± 25.0 events from SM ti pro- 
duction and 13.9 ± 7.5 events from non-tt backgrounds 
in the 2.2 fbr 1 CDF II data sample. 

The mass of the H + candidate is directly reconstructed 
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using the two jets. The mass resolution is improved by 
reconstructing the ti event as a whole with a kinematic 
fitter used for the precision top quark mass measurement 
described in Ref. [13]. The original kinematic fitter is 
modified for the H + search. In the fitter, the lepton, the 
missing Ex (from a neutrino), and the four leading jets 
are assigned to the decay particles from the ti event, and 
the quality of the assignment is evaluated using this \ 2 : 
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The x 2 1S minimized by constraining leptonic W fi- 
nal state (ly) to have the W invariant mass (80.4 
GeV/c 2 ) [16| and both top quark final states {blv and 
bjj) to have the same top quark mass of 175 GeV/c 2 . No 
constraint is imposed on the dijet mass of the hadronic 
boson (jj). In the mass constraints, the transverse ener- 
gies of the final-state objects (p^ mcas ) are allowed to vary 
within measurement uncertainties (oj). The unclustered 
energy (p^ E ' mcas ) is the sum of measured transverse en- 
ergies not included in the leading jets Ex and is used to 
correct the missing Et- In the jet assignment, 6-tagged 
jets are assigned to the 6-quarks. The jets assigned to 
the 6-quarks are called 6-jets, and the other two jets are 
called h-jets. If the tt event has more than two 6-tagged 
jets, the jets with the best \ 2 are assigned to 6-quarks. 
Then, we reconstruct the mass of hadronic boson using 
two h-jets with fit energies (p^ fit ). In this kinematic event 
reconstruction, only 55% of the SM ti events have cor- 
rectly matching jets. The wrong jet-parton assignments 
dominantly come from hard radiation jets which are se- 
lected as leading jets and from the falsely fe-tagged jets 
originating from the hadronic decays of W bosons. 

The expected dijet mass distributions of H + and W in 
top quark decays are produced using the pythia gener- 
ator [ITj and the full CDF II detector simulation. The 
ALPGEN generator [3] with the pythia parton shower 
simulation is used for non-tt backgrounds. In the simu- 
lation sample, the H + is forced to decay solely into cs 
with zero width and with masses ranging from 60 to 150 
GeV/c 2 . 

The simulation shows that the reconstructed H + has a 
significant low-mass tail, which is predominantly caused 
by final-state gluon radiation (FSR) from the hadronic 
decays of the Higgs boson. The hard FSR results in more 
than four final-state jets in a lepton+jets ti event. To 
recover the energy loss due to the FSR, the fifth most 
energetic jet is merged with the closest jet among the 
four leading ones if the pair has a AR distance smaller 
than 1.0, provided that the fifth most energetic jet has 
E T > 12 GeV and |7y| <2.4. Merging the fifth jet results 
in better jet energy resolution and improves the to#+ 



resolution by approximately 5% in more than four final 
jets events for 120 GeV/c 2 Higgs sample. 

In the CDF II data sample of 2.2 fb _1 , we observe 
200 ti candidates in the lepton+jets decay channel. No 
significant excess is observed in the dijet invariant mass 
of top quark decays. Figure [T] shows that the observed 
dijet mass distribution agrees with the SM expectations. 
Hence, we extract upper limits on B(t — > H + b) using a 
binned likelihood fit on the dijet mass distribution. 
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FIG. 1: Observed dijet mass distribution (crosses) compared 
with background distributions of W bosons (filled) and non-tt 
processes (cross hatched) in CDF II data sample of 2.2 fb _1 ; 
the background distributions are added on top of each other. 
An example of the dijet mass distribution from 120 GeV/c 2 
H + bosons (bold line) is overlaid assuming B(t — ► H + b) = 
0.1, which is about the 95% C.L. upper limit on B(t — > H + b). 



The binned likelihood (LH) function is constructed em- 
ploying Poisson probabilities: 



x e 



x G(N hig ,a N ). 



(2) 



The probability of finding events in the mass bin i 
comes from a set of simulated dijet mass distributions 
of H + , W, and non-tt backgrounds. These distribu- 
tions are called templates. The Poisson probability (P l ) 
in each bin is computed from the number of observed 
events, rii, and from the number of expected events, 
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where 



N H + , Nw, and iVbkg are parameters representing the to- 
tal number of events in each template category. The 
minimization of — In C gives the most probable values 
for N H +, Nw, and iVbkg- In the LH fit, N H + and Nw 
are free to vary, however, the non-tt background (iVbkg) 
is estimated independently and is allowed to vary within 
its Gaussian uncertainty (oN bkg ). Based on the number 
of events from the LH fit, a Bit — > H + b) is extracted 
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assuming B(H+ — ► cs) = 1. In Figure [IJ dijet mass dis- 
tributions of the SM events are normalized by the like- 
lihood fit to the observed dijet mass distribution with 
B(t -► H+b) fixed to 0. 

The sources of systematic uncertainty in the extracted 
B(t — > H+b) include uncertainties in the jet energy scale 
corrections flpj j . initial state and final-state radiation, 
modeling of the non-tt background, choice of event gen- 
erators in simulation. These systematic sources perturb 
the shape of the dijet mass and cause a shift in the result 
of the LH fit. The shift in the resulting B(t -> H + b) is 
estimated using "pseudoexperiments" of the perturbed 
and unperturbed dijet mass distributions for each sys- 
tematic source; the pseudoexperiments are generated by 
the bin-to-bin Poisson fluctuations of the simulated di- 
jet mass distributions. The dominant systematic uncer- 
tainty originates from the choice of event generators in 
the simulation, unless rnu+ is close to mw, in which case 
the jet energy scale uncertainty dominates. The other 
systematic uncertainties from data/Monte Carlo differ- 
ences in &-tagging rates and top quark mass constraints 
in tt reconstruction are negligible compared to the un- 
certainties from the perturbed dijet mass shape. 

The individual systematic uncertainties are com- 
bined in quadrature. The total systematic uncertainty 
(AB(t — > H + b)) is represented by a nuisance parameter 
which adds to the branching ratio and has a Gaussian 
prior PDF with width AB(t -> H+b). We eliminate this 
nuisance parameter by Bayesian marginalization (lij and 
obtain a posterior PDF in B(t — > H+b) assuming a uni- 
form prior PDF in < B(t -> H+b) < 1. The expected 
upper limits on B(t — > H+b) with 95% C.L. are derived 
from a thousand pseudoexperiments using the SM back- 
grounds events for each rrifj + . 

The upper limits on B(t -> H+b) at 95% C.L. show 
a good agreement between the observation and the SM 
expectation. The upper limits in Figure [2] includes the 
systematic uncertainty in B(t — > H+b). Since the LH 
fit has very little sensitivity for m H + rs mw, the upper 
limits around 80 GeV/c 2 H+ are omitted in the Figure. 
The exact values of the upper limits in the Figure [2] are 
listed in Table H 



m H+ (GeV/c 2 ) 


60 


70 


90 


100 


110 


120 


130 


140 


150 


Expected 
Observed 


0.13 
0.09 


0.19 
0.12 


0.22 
0.32 


0.15 
0.21 


0.13 
0.15 


0.12 
0.12 


0.10 
0.08 


0.10 
0.10 


0.09 
0.13 



TABLE I: Expected and Observed 95% C.L. upper limits on 
B(t -» H+b) for H+ masses of 60 to 150 GeV/c 2 . 



This analysis can set model-independent limits for 
anomalous scalar charged bosons production in top quark 
decays. Besides the assumption that a scalar boson de- 
cays only to cs with zero width, no model-specific pa- 
rameter is used in this analysis. Therefore any generic 
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FIG. 2: The upper limits on B(t -» H+b) at 95% C.L for 
charged Higgs masses of 60 to 150 GeV/c 2 except a region for 
m H + ~ mw- The observed limits (points) in 2.2 fb _1 CDF 
II data are compared to the expected limits (solid line) with 
68% and 95% uncertainty band. 



charged boson would make a secondary peak in the di- 
jet mass spectrum if it decays into a dijet final state like 
the H+ — > cs in top quark decays. Here, we extend the 
search below the W boson mass (20| down to 60 GeV/c 2 
for any non-SM scalar charged boson produced in top 
quark decays, t — ► X+(-^r ud)b. This process is simu- 
lated for the CDF II detector and is similar to H+ — > cs. 
In the simulation, we obtain a better dijet mass resolu- 
tion for ud decays than for the cs decays. The difference 
in the mass resolution comes from the smaller chance of 
false 6-tagging from light quark final states of X+ than 
the cs decays, thus resulting in a smaller ambiguity of 
jet-parton assignments in the ti reconstruction. Conse- 
quently, the upper limits on B(t — > X+(-^ ud)b) are lower 
than the limits on B(t — » H+(-^ cs)b) regardless of the 
charged boson mass. 

In summary, we have searched for a non-SM scalar 
charged boson, primarily the charged Higgs boson pre- 
dicted in the MSSM, in top quark decays using lep- 
ton+jets tt candidates. This is the first attempt to search 
for H+ — > cs using fully reconstructed charged Higgs 
bosons. In the CDF II data sample of 2.2 fb — 1 , we find 
no evidence of charged Higgs bosons in the dijet mass 
spectrum of the top quark decays. Hence, upper lim- 
its on B(t -> H+b) with 95% C.L. are placed at 0.1 to 
0.3 assuming of B(H+ — * cs) = 1.0 for charged Higgs 
masses of 60 to 150 GeV/c 2 . This analysis also yields 
conservative upper limits on any non-SM scalar charged 
boson X+ production from top quarks. Based on simula- 
tion, we find that the upper limits on the branching ratio 
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B{X + — > ltd) are always better than the upper limits on 
B(JT+ -> cs). 
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